SOME RESULTS FOR POISONING IN A CATALYTIC MODEL 

JEFFREY E. STEIF AND AID AN SUDBURY 
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f^ ' Abstract. We obtain new results concerning poisoning/nonpoisoning in a cat- 

CN , alytic model which has previously been introduced and studied. We show that 

5-H ' poisoning can occur even when the arrival rate of one gas is smaller than the sum 

of the arrival rates of the other gases, and that poisoning does not occur when all 

gases have equal arrival rates. 
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r^ • 1. Introduction 

' Grannan and Swindle (1991) introduced a collection of interacting particle sys- 

tems to model catalytic surfaces. In these models n types of gases represented by 
states {1,2, ...,n} fall on sites of the integer lattice Z'^. A vacant site is represented 
by state 0. Molecules of gas i fall upon vacant sites at rate pi, where '^Pi = 1. No 
(N| ! two different gases can occupy adjacent sites and so if a molecule of type i falls upon 

•^ I a vacant site adjacent to a molecule of type j,j ^ i, then the two gases react and 

— j. ■ both sites are left vacant. If there are several adjacent sites with molecules different 

O . from the arriving molecule, one of them is selected uniformly at random to react with 

the arriving molecule. 

Grannan and Swindle (1991) defined "poisoning" as the configuration converging 
a.s. as time goes to infinity. In such a case, the limit is necessary the "all i configu- 
ration" for some i G {1, 2, ..., n}. Let 6i be the point mass at the configuration where 
everything is in state i and define coexistence as the existence of a stationary distribu- 
tion which is not a convex combination of the {5j}'s. Heuristically, poisoning implies 
^ ■ noncoexistence (but this statement is not precise). Grannan and Swindle (1991) 

proved the following result. (They proved somewhat more than what we state.) 
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Theorem 1. (Grannan and Swindle) (i). Ifn = 2,dis arbitrary and pi ^ 1/2, then 

coexistence does not occur. 

(a). There exists eo > such that for any n and d, if pi > 1 — eo, then for any initial 

state having infinitely many O's, a.s. poisoning occurs with the limiting configuration 

being all I's; i.e., r]t -^ 1, the configuration consisting of all Vs. 

(Hi). For any n and d, if Pi < 2dA L'l+i ("^here \c{d) is the critical value for the 

contact process in d- dimensions), then, for any initial state having infinitely many 

O's, a.s. the system does not converge to the all 1 configuration. 
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(iv). For every d, there exists N^ such that if n > N^, there exist pi, . . . ,p„ so that 
coexistence occurs. For d = 1, one can have coexistence with n = 5. 

(i) was proved using a type of "energy" argument, (ii) was quite involved and (iii) 
and (iv) were proved by arguments involving stochastic comparison with a supercrit- 
ical contact process. It is not known whether coexistence can occur for less than 5 
gases in 1 dimension. 

Mountford and Sudbury (1992) strengthened part (ii) (using a submartingale 
argument) of the above result by showing: 

Theorem 2. (Mountford and Sudbury) Let n and d be arbitrary andpi > 1/2. Then, 
if the initial state,riQ, has infinitely many 1 's or O's, then poisoning occurs with rjt ^ 1 
a.s. 

A question left open is whether it is possible for a gas to poison a surface with 
its arrival rate being less than 1/2. Our first result proves this to be the case. 

Theorem 3. If d = 1 and n = 4, if the first gas has rate 0.47 and the other three 
have rates 0.53/3, then the first gas will poison the surface with probability 1 if the 
initial state of the surface has an infinite number of O's or 1 's. 

It is trivial to observe that if poisoning occurs from every initial state, then 
there cannot be coexistence. However, one should not think of these things as being 
synonymous as we will see. In addition, the proof of Theorem |21 shows that for any d 
and n, if pi > 1/2, then coexistence does not occur. For the interesting case pi = p2, 
Mountford (1992) proved the following result. (What we state is just a special case 
of what he proved.) 

Theorem 4. (Mountford) For d = 1 and n = 2, the only translation invariant 
stationary distributions are {p6i + (1 — p)52}pe[o,i]- 

This strongly suggests that coexistence does not occur here. Our second result 
says that this noncoexistence result is not due to poisoning. 

Theorem 5. Let n and d be arbitrary. If there is i ^ j with pi = pj, then, starting 
from all O's, the probability of poisoning in states i or j is 0. In particular, if all of 
the Pi 's are equal, then, starting from all O's, the probability of poisoning is 0. 

Theorem ini is proved in Section 2, Theorem is proved in Section 3, some further 
comments and questions are given in Section 4 and finally there is an appendix which 
provides some assistance for one of the proofs. 

2. Poisoning in 1-dimension with rate < 1/2 

In this section, we prove Theorem El 

Proof of Theorem 3: The method used will be to consider a block of adjacent I's and 
find conditions under which they tend to spread. At the right-hand end of such a 
block there must be a 0. We condition on the types occupying the next three spaces 
to the right. There are 18 essentially different possibilities as the reader can check; 
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see column 1 of Table 1 for a list of these. Blocks such as 030 and 020 or 203 and 
304 are considered equivalent. 

Consider the same model but defined only on the nonnegative integers. Equiva- 
lently, at time 0, there is a in position —1 and no particles are allowed to arrive at 
this position. We show that if there is a 1 at the origin, then with a fixed positive 
probability, independent of the configuration to the right of the origin, the origin 
remains in state 1 forever and the block of I's containing the origin approaches oo. 
(Of course if the origin remains in state 1 forever, it must be the case that the block 
of I's containing the origin approaches oo.) If this can be done, it follows that in 
the original model, any 1 will spread in both directions poisoning all of the integers 
with a fixed positive probability. Since there are initially infinitely many O's or I's, 
it easily follows that poisoning in state 1 occurs a.s. 

In order to prove the statement concerning the spread of I's, for each configu- 
ration f] on the nonnegative integers which has a 1 or at the origin, we define its 
weight W{ri) as follows. Let B{ri) be the block of I's starting from the origin in the 
configuration r]. (This may be empty or infinite.) If \B{ri)\ = oo (equivalently r] = 1), 
then we let W{ri) := oo. Otherwise, W{ri) is defined to be |-B(^)| plus the score, as 
defined in column 2 of Table 1, of the block following the first after B{r]). The idea 
is to define the scores of the blocks (abc) in such a way that the expected change in 
the weight is always positive, in which case we will obtain a submartingale. 



Table 1 
Scores for essentially different blocks 



block 


Score 


Expression 


Follower 


222 


0.000 


0.0007 


2 


220 


0.163 


0.0012 


2 


202 


0.295 


0.0022 


2 


203 


0.339 


0.0002 


3 


022 


0.354 


0.0034 


2 


200 


0.404 


0.0002 


2 


201 


0.493 


0.0018 


00 


020 


0.498 


0.0031 


2 


002 


0.570 


0.0032 


2 


000 


0.664 


0.0055 


2 


001 


0.827 


0.0058 


00 


010 


0.920 


0.0044 


2 


102 


1.008 


0.0034 


22 


100 


1.157 


0.0036 


22 


Oil 


1.173 


0.0060 


00 


101 


1.456 


0.0056 


00,02 


110 


1.555 


0.0040 


222 


111 


1.997 


0.0054 


0222 



Letting {i]t}t>o denote our process on the nonnegative integers, a very long and 
tedious calculation, left to the reader, shows that there exists c > such that for all 
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rj with B{r]) ^ 0, 

It '*=°-'- ^^> 

[While this long detailed calculation is being left to the reader, the appendix contains 
a discussion which aids the reader in making this calculation; perhaps it takes 2 hours 
of work to check the above with the aid of the appendix and 5 hours otherwise.] It 
follows that with positive probability, uniform in rj, with B{ri) ^ 0, the block of I's 
will grow to cxD before the 1 at the origin changes which is what we wanted to show. 
To carefully do this, we follow the argument for Theorem 3 in Sudbury (1999) and 
proceed as follows. Let T be the (possibly infinite) stopping time when the block of I's 
at the origin disappears. Equation (jT)) implies that {Vr(r7(AT)}t>o is a submartingale 
with respect to the natural filtration of {?7i}t>o- Moreover, it can be shown that 
Equation (^ implies that for e sufficiently small, for any initial configuration ri with 
B{r]) 7^ 0, t/i := 1 — (1 — e)^*^^"^^^ is a bounded submartingale and thus tends a.s. to 
a limit f/oo. Let us assume further for the moment that |-B(r7o)| > 2. This gives us 
that ?7o > 1 - (1 - e)2. If the block of I's dies out, then t/oo < 1 - (1 - e)^-^^^. Since 
i?[t/o] = i?[[/oo], it follows that it cannot be the case that the block of I's dies out 
a.s. Since however Vt converges a.s., it must be the case that the block of I's grows 
to infinity with a uniform positive probability, independent of ?7o with \B(rii^\ > 2. 
If |S(r7o)| = 1, it is clear that with a uniform positive probability, independent of rjo, 
the single 1 at the origin spreads and reaches size at least 2 at time 1. At that point, 
one can apply the previous argument. ■ 

Remarks: 

(i). The score {abc) represents whether this triple is likely to aid the I's in spreading 
or make them more likely to contract. 

(ii). The above proof is not computer assisted in the sense that one can check the 
correctness of the proof without a computer; it suffices to use a hand calculator or 
in fact even hand calculations suffice (the latter requiring a good deal of patience). 
Nonetheless, a computer was essential in helping us find the proof and in particular 
helping us find the scores for the various blocks of length 3. More discussion concern- 
ing this point follows below. 

(iii). One might hope that we could have carried out the above proof using blocks of 
length 2 instead of blocks of length 3 but it seems that this is not possible if we want 
to have pi < 1/2 as also described below. 

While it is not needed at all for checking the correctness of the proof of Theorem El 
we now nonetheless explain in some detail how the scores in Table 1 were arrived at 
which in turn allowed us to obtain the proof. We first tried blocks of size 2. Our 
method was to try various values of pi. As in the proof, we considered only the 
situation at the right-hand end. We wrote down the equations for the instantaneous 
rates of change of the expected weight for the various blocks of size 2. The basic idea 
was to find the values of the (ab), ah ^ 22 which make all these rates 0. Having done 
this we wrote down the equation for (22). If the rate was > then we would have 
found a value of pi which allows the block of I's to spread. If the rate was < 0, then 
we would need a larger value of pi. (In fact we found values of the (ab) which gave a 
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rate very close to and then increased pi so that the rate of change of the expected 
weight was strictly positive for all blocks). When we computed this change in the 
expected weight, we computed it assuming to the right of the block one had what one 
would guess is the most disadvantageous scenario for spreading of the I's. The only 
case considered was P2 = Ps = P'i- A little experimentation with other possibilities 
suggested they would give less favorable results. 

Unfortunately it was impossible to make all the necessary expressions positive 
with pi < 0.5. This meant we had to go up to the next level, looking at blocks of 
length 3 to the right of the block of I's and their immediate 0. There were then 18 
fundamentally different expressions that must be made > 0. Again, we assumed that 
P2 = Ps = Pa- We assumed that 1 > > 2, 3, 4 in the sense that the score for every 
block is increased if the value in the block is increased by replacing it by a "higher" 
value. For example, we assumed that (020) < (000) < (010). This assumption is 
easily checked by looking at the 18 scores in column 2 of Table 1. Assuming this 
worse case scenario to the right of the 3 block, one obtains the following rates of 
change. The assumed worse case scenario to the right of a 3 block is listed in the 
fourth column in Table 1. 
000 

pi((002) + 1 + (100) + (010)) + 3p2(-l + (200) + (020) + (002)/3) + 2^2(000) - 3(000) 

203 

Pi((002) + .5((002) + (200))) +P2(-1 + (020) + (200) + (002)) 
+P2(-1 + (020) + (002) + ((002) + (200))/2) - 2(203) 



100 



111 



001 



Pi{2 + (022) + (110) + (100)) + 3p2((010) - 1 + (000))/2 
+3p2(000) + 2p2(100) + P2(102) - 3(100) 

Pi4 + 3p2(-l + (011)2)/2 - (111) + (1 - pi)((110) - (111)) 



Pi(l+(010) + (101) + (011))+3p2(-l + (000) + (201) + (000))-3(001) + (l-pi)((000)-(001)) 
102 
Pi(2 + (100)) + 3p2(-l + (010) + (002))/2 + P2((002) + (100) + (002)) - 2(102) 

220 

pi((020) + ((200) + (220))/2) +p2(-l + (022)) 
+P2(-1 + (022) + (020) + (220) + (200)) - 2(220) 
022 

pi(l + (002)) + 3p2(-l + (002)) + 2^2(002) - 2(022) 
202 
Pi((002) + ((200) + (002))/2)+p2(-l+(020))+p2(-l+(020) + (002) + (002) + (200))-2(202) 

200 

j>i((000) + (000) + (200)) +p2(-l + (020) + (220) + (202)) 
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+P2(-1 + (020) + (000)3 + 2(200)) - 3(200) 
020 

pi(l + (202) + 1.5(000) + .5(020)) +p2(-l + (002) + (220) + (022)) 
+2p2(-l + (002) + 1.5(000) + .5(020)) - 3(020) 
201 

Pi((001)2)+p2(-l+(020) + (200)-l+(020)+2(001) + (200))-2(201) + (l-pi)((200)-(201)) 
010 

pi{{l + (102) + (110) + (010)) + 3p2(-l + (001) + (000)) +p2((000)2 + (010)) - 3(010) 
002 

Pi(l+(022) + (102) + (000))+p2(-l + (000) + (202) + (022))+2p2(-l + (000) + (302) + (000))-3(002) 
Oil 

pi{l + (110) + (111)) + 3p2(-l + (001) + (001)) - 2(011) + (1 -pi)((010) - (Oil)) 
110 

Pi(3 + (110)) + 3p2(-l + (Oil) + (010))/2 + P2(2(100) + (110)) - 2(110) 
101 
Pi(2+(102) + (lll))+3p2(-l + (010) + (001) + (001) + (100))/2-2(101) + (l-pi)((100)-(101)) 

222 

j9i(022)+j92(-2 + 3(022)) 

A computer program was written which calculated these expressions for blocks 
of any length. The program goes through all possible sets of right-hand blocks made 
up of 0, 1, 2, 3, 4, finds the O's and in turn replaces them with 1, 2, 3, 4. The score of 
each block was then replaced by the score that makes the relevant expression above 
equal 0. For example, the third to last equation above becomes 

(110)„e«, = [Pi(3 + (110)) + 3p2(-l + (Oil) + (010))/2 + P2(2(100) + (110))]/2. 

The first 17 equations were then iterated until the differences between the new 
and old scores given to the blocks were sufficiently close to 0. All 18 equations 
were not used because, since the total rates in and out of the blocks must balance, 
the set of equations would be singular. The last expression, which is for (222) was 
then tested. If it was positive, a suitable set of values would have been found. If 
it was not, pi was replaced by a larger value and the procedure was started again. 
A suitable set of scores was found for pi = 0.4699. The expressions above were 
then independently checked (with the package Minitab) using pi = 0.47 to ensure 
all the expressions came out strictly positive. The first few decimal places of these 
expressions are given in column 3 of Table 1. (It should be noted that the computer 
program tests the transitions for all possible scenarios to the right of the blocks of 
length 3 given above since these scenarios can influence the outcome. The expressions 
given above use simply the best guess as to which of the values 0, 1, 2, 3, 4 will give 
the least advantageous result for the spread of the block of I's. The guesses are in 
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fact correct as can be checked.) The checks made by minitab can then also be carried 
out by hand and/or pocket calculator obtaining the rigorous noncomputer assisted 
proof presented above. 

One can next attempt to prove that even lower values of pi can spread with 4 
types or alternatively study this problem with a different number of types. In either 
case, one has to deal with blocks of length greater than 3 and then the equations 
become too complicated to exhibit or to do calculator or hand computations with. 
The computer program used can only continue up until blocks of length 5 as the 
particular version of Fortran used only allows arrays of size up to 4096. Further, one 
can consider the case when there are infinitely many gases. In this limiting case it 
is assumed that any new arrival will react with an adjacent gas unless it is of type 1 
in which case there will be no reaction if the arriving gas is also type 1. Since this 
situation is more favorable to gas 1, it is possible to use only blocks of length 2 to 
show that in this case, gas 1 can spread with a rate of 0.46. 

Theorem El niay then be improved as follows using the computer program men- 
tioned above. Belief in the following theorem thus relies on belief that the program is 
correct. The program gives the same results as hand calculation for blocks of length 
2 and 3. 

Theorem 6. (Computer assisted proof) In one dimension, poisoning will occur for 
rates of gas 1 greater than those appearing in Table 2 (with all other gases having the 
same rate) if the initial state of the surface has an infinite number of O's or Vs. 

Table 2 
Upper bounds to critical values for poisoning 



Number of Gases 


Rate of gas 1 


Rate of other gases 


2 


>0.5 


<0.5 


3 


0.445(0.40) 


0.277 


4 


0.432(0.38) 


0.189 


oo 


0.46(0.37) 


0.000 



Figures in parentheses are the critical values suggested by simulation. The case of 4 
gases is a worse approximation to the correct value than is 3 gases, as the computer 
program only allows the use of blocks of 5 places in this case. 

3. Nonpoisoning 

In this section, we prove Theorem |31 The argument here is very similar to that 
used in Bramson and Griffeath (1989) and Cox and Klenke (2000). 

Recall that in this result, we are starting with all O's. For each i = 1, . . . ,n, let 
Ai be the event that poisoning occurs and that the final state of all sites is i. Clearly, 
the Aj's are disjoint and their union is the event in question. 

Lemma 7. For each i, P{Ai) G {0, 1}. (This lemma does not require any assumption 
on the rates.) 

Given this lemma, we are done as follows. P{Ai) = P{Aj) by symmetry and by by 
Lemma [71 these are each or 1. Since they are disjoint, they each must be 0. The 
last statement of the result follows immediately. ■ 
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Proof of Lemnia\^ Of course, we just need to consider Ai. The idea is that although 
we have stochastic dynamics, we can encode all the randomness we will need to 
drive the dynamics into random variables associated to each lattice point. To drive 
the dynamics, for each location, we need n independent Poisson processes (for the 
arrivals of the n different types of particles) and we will also need some more random 
variables at each location to be used to decide which particle will be reacted with if 
the arriving particle lands next to more than 2 particles of a different type. 

The details of how to do this are as follows. Let {f/jjjg^d be i.i.d. random variables 
each uniformly distributed on [0, 1]. The entire evolution of our process will simply 
be a function of the {f/j}'s (i.e., it will be determined by the {?7j}'s); no further 
randomness will be needed. On the unit interval [0, 1] with the Borel sets and with 
Lebesgue measure, define random variables {Xt}t>Q, {Yk}k>i and {Zk}k>i such that 
these three processes are independent of each other, the first process is a rate 1 
Poisson process, the Y^s are i.i.d. with P{Yk = j) = Pj for j = 1, . . . , n, and the Z^s 
are i.i.d. with each Zk being a uniform random ordering of the 2d neighbors of the 
origin in Z'^. Note, crucially, that the Xj's, Ifc's and Z^s are functions defined on 
[0,1]. 

Note that we know that on some probability space, we can define random variables 
with the above distribution but the point is that we want to define them on the unit 
interval [0, 1] with the Borel sets. It is well known that this probability space is rich 
enough to be able to define these random variables on it. We of course, as usually 
done in probability theory, do not need to explicitly give what these functions are; it 
is only required that they have the right distribution. 

Now, for each i G Z'^, we can consider the random variables {Xt{Ui)}t>o, 
{Yk{Ui)}k>i and {Zk{Ui)}k>i- These are independent for different i (since the t/j's 
are) and have the same distribution as the {Xt}t>o, {Yk}k>i and {Zk}k>i defined 
above (because Ui is uniform on [0, 1]). Note that Xt{Ui) (as well as these others) are 
random variables since they are a composition of random variables. Note that for this 
to hold, it was important that the underlying probability space for {Xt}t>o, {Yk}k>i 
and {Zk}k>i was the unit interval with the Borel sets and not the unit interval with 
the collection of Lebesgue measurable sets. 

We now use these random variables to drive our dynamics as follows. Fix i G 
Z'^. The arrival times of particles at site i will be taken to be the Poisson process 
{Xt{Ui)}t>o- The type of the fcth arriving particle at site i will be taken to be Yk{Ui). 
If the kth particle arrives and there is at most 1 neighboring particle of a different 
type, we of course know what to do. However, if there is more than 1 such neighboring 
particle, then we look at the random ordering Zk{Ui) of the neighbors of the origin 
and have our arriving particle choose, among those neighbors which have a particle of 
a different type, the neighbor i + u with the highest u value according to the ordering 
Zk{Ui), and then i and this chosen site react and become 0. It is clear that this 
generates our interacting particle system. 

Now the event Ai is of course measurable with respect to {Ui}i^z'' and is also 
trivially translation invariant. Since any i.i.d. process is ergodic (see for example, 
Walters (1975)), it follows that P{Ai) e {0, 1}, as desired. ■ 



POISONING IN A CATALYTIC MODEL 9 

4. Further remarks and conjectures 

We believe the following strengthening of Theorem should be true. 

Conjecture 8. Let n and d he arbitrary. If there are i ^ j such that Pi = pj and all 
the other pk's are no larger, then, starting from all O's, the probability of poisoning is 
0. 

By Theorem the probability that we get poisoned in states i or j is 0. It seems 
natural that it should be even harder to get poisoned in one of the other states with 
a smaller pk- The following monotonicity result seems reasonable and would imply 
Conjecture |H1 Of course, this monotonicity result is not necessary for Conjecture |H1 
and it cannot be ruled out that there are no other types of monotonicity results which 
could be used instead to obtain Conjecture |H1 



nj 



Conjecture 9. Consider two systems with the same d and n with vectors {pi,p2, ■ ■ ■ ,P', 
and {pi,P2,Pn) such that pi < pi and pi > pi for i > 2. Then the process of I's in 
the first process is stochastically smaller than the process of 1 's in the second process. 

We point out however that the simplest naive way of trying to prove Conjecture 
El doesn't work. Consider the case where d = 1, n = 3 and our two systems have 
rates (1/3, 1/3, 1/3) and (1/2, 1/4, 1/4). The obvious way to couple these systems is 
to let particles fall at the same location at the same time with the probability of the 
pair {i,j) falling being given by 

P(2,2) = P(3,3) = 1/4,P(1,1) = 1/3,P(2,1) = P(3,l) = 1/12 

and if both particles have a choice with whom to react with, they choose the same 
one. Now, under this coupling, it can in fact happen that there is a some location 
having a 1 in the first system but not having a 1 in the second system. In the following 
realization, we first have the first process getting a 2 and the second process getting 
a 1. After this, in the next 5 steps, we have the same type particle arrive in the two 
systems and the types, in order of arrival, are 2,3,3,3, and 1. What we then see is 

0000000000,0000000000 

0000200000,0000100000 

0000220000,0000000000 

0000220000,0000300000 

0000020000,0003300000 

0000000000,0003303000 

and finally 

0000010000,0003003000 

or 

0000010000,0003300000. 

5. Appendix 

In this section, we aid the reader in verifying Equation (^ by suggesting how this 
should be done. 
(1). The reader should first check that the 18 expressions in Section 2 correspond to 
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the rate of change of the expected weight in what the reader would think is the most 

disadvantageous scenario to the right of the block. In column 4 of Table 1, these 

supposed most disadvantageous scenarios to the right of the block are listed. The 

case 101 is special in that for certain terms, it is assumed that there is a 02 to the 

right and for other terms, it is assumed that there is a 00 to the right. 

(2). Next, one should plug in the scores for the blocks given in Table 1 into all the 

expressions for the blocks and check that all the numbers one obtains are positive. 

The first few decimals for these numbers are listed in column 3 of Table 1. 

(3). Last, one has to consider all the other possible scenarios which can in fact sit 

to the right of the block, compute how each one affects the rate of change of the 

expected weight and check that all these other cases yield a larger value. 

One should observe that our assumed worse case scenario sometimes holds uni- 
formly and sometimes not in the following sense. A 00 to the right of the block 001 
is always worse no matter what arrives next. Similarly, a 2 to the right of the block 
022 is always worse than a 0. However, for the block 000, it is not uniformly worse to 
have a 2 to the right rather than a 0. Should a 3 arrive at the right most 0, it would 
have been worse to have a to the right since then the 3 would have remained but 
with a 2 to the right, the 3 would not stay. However, if a 1 would have arrived at that 
position, it would be worse to have a 2 since then the 1 would not stay. However, in 
all cases, when one averages over all the possibilities, these assumed worse cases are 
in fact worse case. However, this needs to be checked and this is precisely step (3) 
above. 
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